Rotating molecules in optical lattices, alignment and monopole crystals 
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The recent progress towards production of near-ground state quantum-degenerate molecules raises 
the issue of how such "small" molecules behave in an optical lattice. In this Letter we show that 
the coupling of the molecular orientation to the local electric field direction will provide several 
new phenomena. In the case where the lasers forming different crystallographic directions of the 
lattice are incoherent, the orientation of the molecules is conserved (for L — 1) and a novel form 
of anisotropic superfluidity can be expected. When the lasers are coherent, and the optical lattice 
is such that the splitting of the rotational levels is large compared to the centre of mass energies, 
an adiabatic description of the molecular orientation is appropriate. This leads to geometric vector 
potentials, pseudo-magnetic monopoles and a frustrated band structure with degenerate minima. 

PACS numbers: 03.75.Mn, 03.75.Nt 
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The rapid progress in making molecules in ultracold 
atomic gases using Feschbach resonances [l| and more 
recently using photoassociationQ, Q has provided a new 
direction for degenerate gas physics, both in terms of 
non-equilibrium effects due to Feschbach sweeps of mag- 
netic fields but also equilibrium many body phenomena 
such as the BEC-BCS crossover. The first non-s wave 
molecules have been studied in CS2QI and very recently 
in K^jH]. Novel phenomena such as splitting of the dif- 
ferent multiplet Fcshbach resonances has been observed. 

It was soon realised that optical lattices could pro- 
vide a very sheltered environment for the production of 
molecules due to control of occupancy (leading to a lack 
of harmful three-body collisions). Initially this led to 
photoassociative state-selective production|6[. More re- 
cently the combination of the tunability of the optical lat- 
tice in terms of tunneling (thorough laser intensity) and 
tunability of the interaction through Feshbach resonances 
has been used to make molecules with both fermionsQ 
and bosons [1]. The reliability of the occupancy in the 
lattice of the molecules is now very goodQ. 

Very recently Raman down steps (STIRAP)[Io| to- 
wards ground state molecules have been achieved and 
with the proposal of frequency-comb methods^, it is 
timely to consider the nature of near-ground state molec- 
ular quantum states in optical lattices. 

Large molecules (i.e. with an extent large compared 
to the optical lattice parameter) have been studied the- 
oretically in optical lattices by Fedichev et alfllj|. The 
nature of binding on scales large compared to the optical 
lattice parameter were established in some detail. For 
such a large s-wave molecule there is no coupling to the 
orientation of the lattice except through the anisotropy 
of the center of mass motion. 

In this Letter we consider small molecules where there 
is a significant orientational coupling to the optical lat- 
tice due to their inherently anisotropic polarisability. We 
will see that this coupling has several pronounced effects 
on the cenrte of mass motion, raising issues of the nature 
of resulting condensed srtates. (For larger molecules, this 



presumably corresponds to interference between the in- 
cident optical field and light scattered by one of the con- 
stituent atoms, as experienced by the other atom.) 

The most striking effects are in an excited rotational 
state which might cause disquiet in terms of lifetimes of 
such states. However we take comfort from the case of 



ortho- and para-Hydrogen LJ,[15[, where ortho- Hydrogen 
(which has L = 1) has a lifetime at condensed matter den- 
sities of the order of a day due to the very small matrix 
elements for relaxation of the nuclear I tot = 1 triplet to 
the singlet. 

In fields of the strengthpjl, EH occurring in optical 
lattices, L 2 is a good constant of the motion despite the 
manifest lack of spherical symmetry. This is due to the 
rotational level spacing being typically six orders of mag- 
nitude larger than the recoil energy. However multiplet 
splitting (as against mixing) already indicates coupling of 
the orientation to the local field direction and it is the 
consequences of such a coupling in an optical lattice that 
we will explore in this Letter. 

Let the molecule have polarisabilitics a" and a~ L paral- 
lel and perpendicular to the molecular axis respectively, 
and define the average polarisability, a = |(er' + 2a^), 
and 8a = cv" — a 1 -. We assume from now on that the 
lattice is red-detuned and that 8a > 0. Then to second 
order in the applied field, £(z), the contribution to the 
Hamiltonian, 7i"\ for a given multiplet is: 
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(1) 

(see, for example ref ([16[) where the angular momentum 
quantum numbers are I and m with the axis of quantisa- 
tion parallel to the field. 

For a one-dimensional position-dependent electric 
field, £(X)z, where X is the centre of mass position of 
the molecule, we define (now considering L = 1): 

H Z (X) = -\£ 2 {X){a 1 +a 2 {m)\m z =0)(m z = 0|) (2) 

where z denotes the direction of the field, ai = a — 
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j^5a and a 2 = ^Sa. The physical interpretation of this 
is semiclassical. For the states m = ±1 the molecule 
presents its less polarisable directions to the field. This 
means the barriers to tunneling arc reduced compared 
to random orientation but also the minima are raised 
compared to random orientation, m = presents (semi- 
classically) the most polarisable axis half the time and 
hence has higher barriers, but lower minima. 

We see that in one dimension the orientation of the 
molecule is a constant of the motion. It is perhaps more 
surprising that we can make statements of a similar na- 
ture in higher dimensions. However the mutual coherence 
or incoherence of the lasers in the different axes of the 
optical lattice affects the results substantially. 

In the incoherent case, let the polarisation of the elec- 
tric field of the laser associated with optical lattice axes 
i = x, y or z be denoted by the unit vectors ej. Then the 
perturbed energy is the sum of three terms of the form 
©: 

W(R) =H & » (X) + He y (Y) + % s (Z) 

=-\ Yl £2 ( R ^ ( ai + a ^ m - = °>( m * = °i)( 3 ) 

i— x,y,z 

where m, is an eigenvalue of • L. 

The result is particularly straightforward if the set {e^} 
is orthogonal. Then the set {\rrii = 0)} (the "polar ba- 
sis" ) is also orthogonal. In that case the three contribu- 
tions to the Hamiltonian from the projection operators 
\rrii = 0)(mj = 0| commute and hence the Hamiltonian is 
diagonal. A molecule in one of the eigenstates, |m, = 0), 
will experience two directions with small periodic poten- 
tials and in the direction parallel to e^, a deeper one. For 
L > 1, the situation is more complicated: for example it 
is not guaranteed that the m = states along three axes 
are orthogonal (e.g. L = 2), so the projection operators 
do not commute. 

The basis set allows some translation of the S = 1 
spinor results from references [3] . The molecule-molecule 
scattering will have singlet and quintuplet scattering 
lengths which determine whether the ground state is 
a polar state (as in the 5 = 1 case) which is com- 
patible with the diagonaliscd Hamiltonian, albeit with 
anisotropic supcrfluid properties. It is less clear what 
happens for the case which would be ferromagnetic with- 
out the lattice. 

Let us now turn to the coherent case. Here the polar- 
isation of the optical lattice is position-dependent. For 
definitcness, pick 

E(R) = £ (x sin kZ + y sin kX + z sin kY) cos Qt (4) 

and for future reference define |E(R)| = £(R). Then the 
antinodes, i.e. the minima of the optical lattice potential, 
are at R = (A/4) (2^ + 1, 2m + 1, 2n + 1). There arc four 
distinct minima (as the Hamiltonian docs not distinguish 



the sign of the electric field), due to the relative signs of 
the components of E. 

The case where there is an adiabatic separation of the 
energy scales of splitting the angular momentum multi- 
plet and the centre of mass motion in one of the optical 
lattice wells (and hence in an even more pronounced form 
in terms of tunneling) will be the focus for the remaining 
part of this Letter. In that limit the molecule preserves 
its orientational state relative to the local electric field 
direction, determined by the centre of mass position. Un- 
like the incoherent case the results are qualitatively simi- 
lar for all L > 0. Berry and Robbins[I3] studied a related 
problem where there was a linear coupling between angu- 
lar momentum and spatial variables in an adiabatic limit 
which we will comment on presently 

We form the adiabatic basis, |m,R), by rotating the 
reference basis set, |m), so that the z-axis is transformed 
to the local electric field direction appropriate to centre 
of mass position R: 

|m,R) = D(-f7(R),-C(R),0)|m) 

= cxp(i?7(R)cr z ) exp(iC(R)crj,)|m) 

where {ci}, with i — x, y or z, are the spin-1 Pauli 
matrices. Thus 

E(R) • L\m, R) = m\m, R) 

with E(R) is the unit vector corresponding to the local 
electric field. 

We then look for an approximation to the solution, 
^(R, C, ?7), of the Schrodinger equation for the molecule 



^V 2 R - i£ 2 (R) (ai + a 2 |0,R)(0,R|) 
The adiabatic approximation is to assume 



*(r,0?)~ i(a)KR) 



which leads (upon substitution and performing an in- 
tegration over the angular variables) to the effective 
Schrodinger equation in terms of a vector potential ppj 
A nl = -i(n,R|V R U,R) : 



A' li • A lm + V nm (R.) ) 4> m CR) = Eif; n (R) (5) 

2m 



where wc define the covariant derivative, A. nl = S Vr — 
iA nl and 

V nm (R) = -j£ 2 (R) {ax5 nm + a 2 Po m )+J2 Anl - Almsnm 

(with no summation convention), where P^ 171 = 6 m0 6 n0 . 
Non-Abclian vector potentials lead to effects experienced 
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by multi-level atoms in light fields: the case of spin- 
orbit effects plj. laser- assisted tunneling and tripod- 
field induced degenerate dark states [23}. This Letter con- 
tains the first example the new phenomena experienced 
in a t/iree-dimensional extended setting with such non- 
Abelian fields. 

We may evaluate the vector potential: 

A"' = [{af cos C + of sin C)Vrt? + VrC] 

and hence determine the associated field: 



l [A-i, Aj 



= 2(d i r l d j ( - 9^00™ cosC - cr™ sinC) 

implying that the pseudo -field, B" m = is 

B mn = 2VC x Vr](a™ n cos ( - a™ n sin (6) 

with the polar angles, £ and 77, playing the role of "Eulcr 
potentials" for the field B. 

An immediate consequence of equation ^ is that V • 
B = except possibly at singularities of V77 and/or V£, 
which occur at the nodes of the electric field of the optical 
lattice (@|. Firstly consider the form of the electric field 
in the vicinity of the node at the origin 

E(R) ~ kS {xZ + yX + zY) 

To examine the nature of the pseudomagnetic field and 
its flux, it is easiest to calculate ^ in polar coordinates 
around the node. The relation between the polar coordi- 
nates of the Electric field vector, (£, 77), and those of the 
centre of mass {9, <p) are determined implicitly from: 

E x ~ sin C cos 77 = Z/R = cos 9 

E y = sin C sin 77 = Xj R = sin 9 cos <fi 

E z = cosC = Y/R = sin 6* sinc/> 

Using expressions for cos £ and tan 77 we find after some 
algebra that: 



VC x V77 



R 1 

^sm~C 



and hence 



W nn (R) 



9 R 



cot C) 



(7) 



in the vicinity of the node. 

We will now calculate the flux through a spherical sur- 
face of B m ™. We will neglect the Dirac-likc strings which 
emerge along the positive and negative y-ax.es (reminis- 
cent of the two strings in Schwinger's treatment of 
monopoles), as the strings are irrelevant to our final lat- 
tice treatment. 

The second term in Eqn. ([7]) vanishes as there is a 
cancelation from the contributions from < 6 < it and 



7r < <fi < 2ir. The first term provides 8wa z flux through 
the surface, corresponding to monopoles of charge ±87r. 
These are the analogues of the monopole in the the work 
of Berry and Robbins[l7j. 

The field in the vicinity of the nodes at R = (m, n, £) 
has the form: 

E(R) oc (-l/(Z-^)x+(-l) m (A-m)y+(-l)"(r-7i)z 

In the case of one or three of £, m and 77, being odd, then 
either the sign of V77 or VC reversed. Thus the (m = 1) 
monopole has negative charge. Conversely, if there are 
zero or two odd integers, the (m = 1) monopole has a 
positive charge as neither or both gradients are reversed. 

Thus the monopoles reside in a Na-Cl lattice, coinci- 
dent with the nodal points of the electric field, of alter- 
nating positive and negative magnetic charges. The signs 
of all charges are reversed if the angular momentum of 
the molecule is reversed. 

Since the adiabatic limit will require a strong optical 
lattice it is natural to construct the tight binding model 
corresponding to the continuum Schrodinger equation, 
Eqn. ([5]). The lattice is cubic and in the adiabatic limit 
the distinction between the field directions on the dif- 
ferent sites in the unit cell is irrelevant as they are all 
perfectly adiabatically connected and degenerate in en- 
ergy. The monopoles reside on the dual lattice, with the 
opposite signs adopting a Na-Cl structure. We will focus 
on the m — 1 case. 

We must assign Peierls factors, e 1 " 4 "."' , to the hopping 
term on each link between sites n and n'. The flux, 
$ 



n' .n" .n" 

* _ 



as e 



/ , through a plaquette is counted modulo 2tt 
Then we must pick A n n > so that we 
get a flux of ±27r/3 through each plaquette, so we ob- 
tain a total flux out of the six plaquettes surrounding 
a monopole of =F87r. The Peierls factors are chosen as 
indicated in Fig. (p}, where we see there are two sites 
per unit cell denoted A and B. This leads to the secular 
determinant (here t > and lo = e l27r / 3 ): 



-2t:(cj- 1 A + Y + loZ) 



-2t{uX + Y + uj- 1 Z) 



(8) 

where: X = cos k x and similarly for Y and Z . This 
implies 

c±(k) - ±V2~W(X - y)2 + (Y - Zf + (Z- X f (9) 

The two bands touch along the (±1,±1,±1) directions, 
with the dispersion relation in that vicinity (defining 
fell parallel to the (1,1,1) direction and Ar 1 (> 0) radi- 
ally perpendicular to the direction) being e±(k\k ± ) ~ 
±V3\sin(0/V3)\k- L for fc^/fcH < 1. 

The most important aspect of the dispersion relations 
relates to the frustrated nature of the Hamiltonian. The 
minima of e_(k) are at the three points in the Brillouin 
Zone: k = (0, 0, it) and cubic equivalents. Because the 
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FIG. 1: The choice of Peierls factors in Eqn. ([8)l. An arrow 
from B to A, for example, indicates a factor of e 12 "^ 3 in the 
term hopping from B to A. 



(direct) lattice is FCC, the other candidates k = (0, 7T, ir) 
are outside the first Brillouin Zone and equivalent to the 
first set by adding a reciprocal lattice vector, 7r(l, 1, 1), 
of the BCC reciprocal lattice. Thus the minima in en- 
ergy are —At which is higher than the — 6i expected of 
a simple cubic system. The wavefunctions on the A 
and B sublattices, (ipA,ipB), for the three minima are: 
for k<*> = (0,0, tt), {^\^b ] ) = (I,""); for = 

(o,tt,o), (vi y) ,^T) = for = (t,0,0), 

The three minima, in conjunction with the possible 
competition of singlet and quintuplet scattering lengths, 
make the nature of the mean field ground state in general 
unlcear. 

In conclusion we have shown that in both incoherent 
and coherent optical lattices the orientation of a molecule 
with L / couples to the direction of the optical field. 
In the incoherent case the orientation is conserved in a 
highly symmetric manner for L = 1. For the coherent 
case the strong optical fields lead to the angular momen- 
tum around the local electric field being an adiabatic 
invariant. The consequences of the adiabatic approxi- 
mation failing at the nodal points of the lattice, in the 
classically disallowed region, affects the tunneling of the 
molecule. This leads to pseudomagnetic monopoles with 
their flux influencing the band structure of the molecule 
in a striking manner. The consequences for condensa- 
tion remain to be explored, as does the configuration of 
monopoles in cases of general (eg non-bipartite) dual latt- 
tices. 
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